Interference pattern of a long diffusive Josephson junction 
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We calculate the modulation by a magnetic field of the critical current of a long disordered 
Josephson junction in the diffusive limit, i.e. when the dimensions of the junction are larger that 
the elastic mean free path, and when the length L is much larger than the width w. Due to the 
averaging of the gauge invariant phase factor over diffusive trajectories, the well-known oscillations 
of the Fraunhofer pattern are smoothed out and replaced by an exponential decay at large field. 
The predicted pattern is universal, i.e., it is independent of the disorder strength. We point out an 
interesting relation with the physics of speckle correlations in optics of turbid media. 
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Introduction - The supercurrent flowing through a 
tunnel junction between two superconductors is given by 
the well-known gauge invariant Josephson relation 
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long diffusive junction, the critical current varies as 
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where (j) is the flux through the junction and = h/2e 
is the superconducting flux quantum. 



where 8 is the phase difference between the two supercon- 
ductors. The Josephson effect is thus a beautiful tool to 
exhibit interference effects, manifestations of the phase 
coherence of the superconducting wavefunction. For ex- 
ample, a circuit with two Josephson junctions is a realiza- 
tion of Young's two slits experiment, where the interfer- 
ence is modulated by the Aharonov-Bohm flux through 
the circuit [1]. Moreover a single Josephson junction with 
a finite width exhibits an interference pattern reminis- 
cent of the diffraction (Fraunhofer) pattern of a slit, as 
recalled in eq. (4) [2]. 

It is natural to wonder whether such an interference ex- 
periment can probe phase coherence in a more complex 
medium with multiple scattering of the electrons. Here, 
we consider a long Josephson junction made of a dif- 
fusive metal forming a quasi-one-dimensional wire. The 
junction of length L is attached to superconducting leads 
along the direction x. A magnetic field is applied along 
the direction z perpendicular to the wire. The width of 
the junction (along the y direction) is denoted by w, and 
its width (along z) is denoted by h. We consider a long 
junction such that L ^ w. The junction is schematically 
represented on figure 1. The amplitude of the Josephson 
current in such a diffusive junction has been calculated 
with the Usadel equation [3, 4]. It has been found that, 
contrary to the the case of the tunnel junction, for a good 
contact between the metallic region and the superconduc- 
tors, the Josephson relation may not be sinusoidal. How- 
ever, harmonics are expected to decay rapidly, roughly as 
(— l)"/n^. In this letter, we assume a sinusoidal Joseph- 
son relation and consider how the phase is modified by 
the application of a magnetic field. We find that for a 
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FIG. 1: Geometry discussed in the text. The junction (grey) 
is a diffusive metal so that the current through the junction 
results form the contribution of many diffusive trajectories. 
We assume that the junction has a quasi-lD geometry: L 3> 



Josephson current and diffusive trajectories - Quite 
generally, the Josephson current resulting from all cur- 
rent paths has the form 



/(<5)=/o^sin (^S{C)-^J^A- 



dl 



where (■ ■ ■ )c denotes the average over all current paths 
through the junction. 

We choose a gauge where the vector potential A is 
aligned along the direction x, = By, y G [—w/2, w/2]. 
On figure 1, the diffuse paths 3—4 represent current paths 
C while the straight path 1 — 2 serves as a reference path 
and corresponds to y ~ 0. The circulation of A is zero 
along the paths 1 — 2, 1—3 and 2 — 4, so that the phase 
difference does not depend on y and is denoted Sq. The 
current can be rewritten in the gauge dependent form 
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I{S) - lo {^sin (^So-^ J^A-dl 
so that wc write the critical current Ic = max[I{5)] as 
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Ic = Io 

or, in a gauge independent form: 



4>{C) is the flux through the area 1 — 3 — 4—2 defined by 
a diffusive path C. 4>o = h/2e is the superconducting flux 
quantum. We have neglected the current flowing in the 
superconductors along the 1 — 3 and 2 — 4 segments, i.e. 
we have assumed that the penetration length A of the 
magnetic field in the superconductor A ^ 0. Taking into 
account a finite A would amount to replace L by L + 2A, 
as usual. 

As a reminder, we first briefly consider the case of the 
short ballistic junction. The current has to be summed 
on the current paths 
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leading to the well-known Fraunhofer-like result [2] 
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where cf) = BwL (or (p = w{L + 2X) to account for a finite 
penetration length). 

For a long diffusive junction, the phase factor in (3) 
has to be averaged on the distribution of diffusive tra- 
jectories. In order to perform this average, we need to 
describe the diffusion from a point r at one end of the 
diffusive sample to another point r' located at the other 
end. We introduce the probability P{r, r', t), solution of 
the covariant equation 
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where the electron charge is denoted — e. This solution 
may be expressed as a functional integral [6] : 
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We consider a long junction where the dephasing be- 
tween r and r' is supposed to be independent of the 
position of r and r' on the boundaries, that is indepen- 
dent of the coordinates y and z. Therefore, we consider a 



one-dimensional diffusion equation with the appropriate 
gauge: 
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P{x,x',t) = 6{x-x')S{t) 



from which we obtain the average 
jectories 
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where Pq is solution of eq. (7) with _B = 0, and x, x' 
are taken at the extremities of the junction. Then the 
critical current is obtained from eq. (3). 

We solve this equation with the magnetic field as a 
perturbation. The eigenvalues of this diffusion equation 
are solutions of 
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and are given by 
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The new magnetic field dependent term implies an 
exponential decay of the probability P(x,x\t) = 
Po{x,x' ,t)e~*^^'^ , with the characteristic time tb given 
by [7]: 
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The average (8) over diffusive trajectories is thus related 
to the Laplace transform P^{r,r') = / Po{r,r' ,t)e~'^*dt 
of the probability to diffuse from one end to the sam- 
ple to the other. The numerator is the solution of the 
differential equation: 



(12) 



with 7 = 1/tb, and with the appropriate boundary con- 
ditions. Here we assume that the disordered junction is 
connected to reservoirs and express that the probabil- 
ity vanishes at the edge of the diffusive metal (different 
boundary conditions could be discussed, but lead also to 
the same 1/ svahL/Ls of eq. (15)). The solution of this 
equation is [5] 

„, ,^ _ Lb sinhx,n/£Bsinh(£ - xm)/Lb . 
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where Xm = niin(a;,a;') and Xm = max(a;, a;'). Wc have 
introduced the characteristic length: 



Lb = \/Dtb 
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The coordinates x and x' are close to the end of the 
diffusive junction. Their value is respectively I and L — I 
where the length / is of the order of the elastic mean free 
path If, [8]. Since I <^ L, we obtain 



P^{l,L- I) 
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In the limit Lb oo, the probability scales as l/L which 
expresses Ohm's law that the transmission coefficient of a 
diffusive system scales like the inverse of its length. Now, 
from (14) and (8), we obtain finally a result independent 
of 



L = Io 



L/Li 



sinh L/ Lb 



that we write in the final form (2) where i 
flux through the junction. 



(15) 
BwL is the 
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FIG. 2: Comparison between the Fraunhofer pattern of a 
short ballistic junction (dotted line, eq. 4) and the pattern 
of a long (L w) diffusive junction (continuous line, eq. 2), 
for the magnetic field dependence of the critical Josephson 
current. (f> is the flux through the junction. 



Phase coherence - This result assumes full phase co- 
herence in the metallic junction. Wc now introduce a 
finite coherence time and the probability P^{x^x') is 
now solution of eq. (12) with 7 = \/tb + I/t,/,. We 
immediately obtain 
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where 1/i^ = l/i^ + l/i^ and L^ = \JL)t^ is the 
phase coherence length. Similarly, the effect of a finite 
temperature can be taken into account by the thermal 
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\J~DJT, so that 1 /L\ is replaced by 1 /L\ - 



Gaussian accumulation of the phase - We now try to 
give a simple interpretation of our result. The dephasing 
accumulated along diffusive trajectories is characterized 
by the average (e~''^)c of the phase factor ip = ^ J A-dl 
along all diffusive paths C in the junction. Since diffusion 
is a Gaussian process, the average over trajectories of a 
given length, that is of a given diffusion time t, is 



For a quasi-one-dimensional diffusion, the average ((Pt)c 
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is simply given by (ipf 

average taken along the transverse direction. Since = 
B^w^ /12 and {Xf)c = 2Dt, we immediately obtain that 
the phase factor averaged along all trajectories of time t 
is (e~"^')c = e"*/"^^ where tb has been defined in (11). 
Then the dephasing has to be averaged over all times t 
for trajectories crossing the sample {x ^ l,x' ^ L — I) 
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which is nothing but eq. (8). 



Relation with weak localization - The magnetic field 
dependence of the Josephson current probes the phase 
accumulated along diffusive trajectoires which cross the 
sample. This physics bears of course some similarity with 
the weak localization correction which probes the distri- 
bution of dephasing along closed trajectories. Instead of 
probing the probability to cross the sample P,{1-, L — I), 
the weak localization correction probes the return prob- 
ability P,{x, x). As a result, for large L ^ Lb, the weak 
localization correction to the dimcnsionlcss conductance 
(in units of 2e? /h) decays as l/L. 

Ag = -/ P,{x,x)dx ^ (18) 

L Jo l:$>Lb L 

while in the same limit L ^ Lb, the Josephson current 
decays exponentially 
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This behavior is very reminiscent of the structure of the 
harmonics of the weak localization correction on a ring as 
g-mL/Le f^.^^^ so-called Alsthuler-Aronov-Spivak oscilla- 
tions [9]). It is a signature of the Gaussian decay of the 
probability to diffuse from one end to another after a time 
t, which scales as /^Dt^ That is why the weak local- 
ization correction (0*'*-harmonics) is a power law while 
the TO 7^ harmonics decay exponentially. In the case 
of the ring, the boundary conditions are periodic along 
the ring, leading to the g"™^/^^ decay of the harmonics. 
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Here, the trajectories can diffuse m times back-and-forth 
before leaving the sample, leading to contributions of the 
form e-(2™+i)i/-f'B. The l/sinhi/Ls behavior results 
obviously from the additive contributions of these diffu- 
sive trajectories : l/sinhL/Ls = J2m e-^^^+i)^/^^ . 

Relation with experiments in optics - Diffusing Wave 
Spectroscopy - The Josephson relation (1) involves the 
transmission coefficient of Cooper pairs which carry ran- 
dom phase factors that have to be averaged over diffusive 
trajectories. It is interesting to notice a similarity be- 
tween our result for the diffusive Josephson junction and 
some results obtained in the physics of speckle correla- 
tions in optics. In optics, the so-called Diffusing Wave 
Spectroscopy (DWS) is a technique consisting in mea- 
suring the correlation function of the transmission am- 
plitude t of light through a turbid medium, measured at 
different times and T [10] . If the scatterers of the dif- 
fusive medium can move, the correlation function gives 
some information on the motion on the dynamics of the 
scatterers. 

The product {t{0)t*{T)) involves pairings of diffusive 
trajectories which carry slightly different phases, since 
the scatterers have moved. Therefore it measures an av- 
erage phase factor: 

(t(o)r(r)) = (|t(o)|2) (e-'^)c . 

The phase accumulated along diffusive trajectories de- 
pends on the dynamics of the scatterers. For example, 
for a Brownian motion of the scatterers, the phase fac- 
tor accumulated along trajectories of time t decays ex- 
ponentially, (e~*'''')c = 6"*/"^^, where the characteristic 
dcphasing time depends on the ratio between the wave 
length A of the incident light beam and the typical dis- 
placement y/D^T of the scatterers after time T. It has the 
form ~ T^X^ / DhT, where Tg its elastic mean path and 
Db is the diffusion coefficient for the Brownian motion of 
the scatterers. Consequently the phase factor averaged 
over all trajectories which cross the sample is obtained 
from (17) and has the same decay (16) as found here 
for the diffusive Josephson junction in a field. Angular 
and frequency speckle correlations (the so-called Ci cor- 
relations) exhibit a similar behavior [5, 11]. The common 
physical origin is the Gaussian accumulation of the phase 
along diffusive trajectories. 

Conclusion - comparison with experiments - Surpris- 
ingly, experiments on diffusive long diffusive SNS junc- 
tions are pretty recent, the difficulty being of keeping 
phase coherence along the junction [12]. A recent experi- 
ment have measured the interference between two metal- 
lic Au long junctions sandwiched in a superconducting 
Al circuit [13]. The total current oscillates with the flux 
through the circuit which modulates the relative phase 
between the junctions. This interference pattern is mod- 
ulated by the interference pattern of each junction. This 



modulation is well described by our result. The low field 
behavior is very well fitted without adjustable parameter 
by the expansion Ic = Iq{1 — j^^) + ■ ■ ■ of our result 
(2). At large field however, the decay seems to be faster 
than exponential. 

During the completion of this work, we have been 
aware of a preprint by Hammer et al. who consider the 
critical current of a long diffusive junction, within the 
Usadel formalism. For the case of a perfect transmission 
at the NS interface, they find for large L/ Lb an exponen- 
tial decay of the form L/LBe~^^^'^ which is compatible 
with our equation (15) but they have considered only nu- 
merically the full range of magnetic field [14]. Another 
paper [15] solves Usadel equation in the limit L ^ Lt- 
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